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Abstract. We present here the discovery and characterisation of a very light planet around HD4308. The planet orbits its star 
in 15.56 days. The circular radial-velocity variation presents a tiny semi-amplitude of 4.1 ms"' that corresponds to a planetary 
minimum mass m2sini= 14. IMe (Earth masses). The planet was unveiled by high-precision radial-velocity measurements 
obtained with the HARPS spectrograph on the ESO 3.6-m telescope. The radial-velocity residuals around the Keplerian solution 
are 1.3 ms"', demonstrating the very high quality of the HARPS measurements. Activity and bisector indicators exclude any 
significant perturbations of stellar intrinsic origin, which supports the planetary interpretation. Contrary to most planet-host 
stars, HD4308 has a marked sub-solar metallicity ([Fe/H] = -0.31), raising the possibility that very light planet occurrence 
might show a different coupling with the parent star's metallicity than do giant gaseous extra-solar planets. Together with 
Neptune-mass planets close to their parent stars, the new planet occupies a position in the mass-separation parameter space 
that is constraining for planet-formation and evolution theories. The question of whether they can be considered as residuals 
of evaporated gaseous giant planets, ice giants, or super-earth planets is discussed in the context of the latest core-accretion 
models. 
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^ 1. Introduction 



After a decade of enthusiastic discoveries in the field of extra- 
solar gaseous giant planets, mainly coming from large high- 
precision radial-velocity surveys of solar-type stars, the quest 
for other worlds has now passed a new barrier with the re- 
cent detection s of several p lanets iii the Neptune-mass regime 
jMcArthur e"alj^2004 .Santos et alJ|2004b: Butler et al...20o4 

23 



IVogt et allbOOl iRivera et al Jl2005f iBonfils et al JbOO*! . They 
are supposedly mainly composed of icy/rocky material, being 
formed without or having lost the extended gaseous atmosphere 
expected to grow during the planet migration towards the cen- 
tre of the system. 

This new step forward has been made possible primar- 
ily thanks to the development of a new generation of instru- 
ments capable of radial-velocity measurements of unprece- 
dented quality. The "fer de lance" among them is undoubt- 



Send offprint requests to: S. Udry 

* Based on observations made with the HARPS instrument on the 
ESO 3.6 m telescope at La Silla Observatory under the GTO pro- 
gramme ID 072.C-0488. 



edly the ESO high-resolution HARPS fiber-fed echelle spectro- 
graph especially designed for planet-search programmes and 
asteroseismology. HARPS has already proved to be the most 
precise spectro-velocimeter to date, as it reaches an instru- 
mental radial-velocity accuracy a t the leve l of 1 ms"' over 
months/years (Mayor_et_aL.20Q3.: Lovis et alA2005). The preci- 
sion a chieved is even better on a short-term basis ( iBouchv et'al] 
|2005'). Another fundamental change that has allowed this 
progress in planet detection towards the very low masses is the 
application of a careful observing strategy to reduce as much 
as possible the perturbing effect of stellar oscillations hiding 
the tiny radial-velocity signal induced on solar-type stars by 
Neptune-mass planets. This is discussed further in detail in 
Sect.E 

In about 1.5 years, the GTO HARPS planet-search pro- 
gramme has discovered (or co- discovered) 1 1 planets, amongst 
which ther e are: 5 h ot Jupiters JPepe et al 2004-'Mou tou et alJ 
2005: L o Curto et al.. 2005) : 3 planets with masses in the sub- 
Saturn mass regime and orbiting at m oderate distances from 
their parent stars ( Lovi s et al and 3 Nept une-mass 

planets, namely a 14.4 Mg planet around yuAra dSantos et alJ 
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Fig. 1. Distribution of planet to primary-star mass ratios for the 
known exoplanets and low-mass brown dwarfs. The HARPS 
detections are represented by the hatched histogram. Planet 
symbols indicate the positions of Jupiter-, Neptune-, and 
Uranus-mass objects. 



| 2004bl) . a 16.6Me planet around the M3 G1581 dBonfils et alJ 
|2005 j. and the planet described in this paper As illustrated 
in Fig^ these planets lie mainly in the low-mass tail of 
the planetary-mass distribution. This distribution, which was 
known to be biased towards low masses, can now be explored 
in greater detail. The characterisation of the low-mass objects 
will strongly constrain planet-formation and evolution scenar- 
ios, as planets with masses between 10 and 100 Me are not ex- 
pected in large numbers according to some current formation 
models (see e.g. Ilda & Linl2004al) . 

The discovery of the extremely low-mass planets represents 
a new benchmark for planet surveys and provides information 
on the low end of the planetary-mass distribution. In this paper, 
we present the discovery of a new very low-mass companion 
to the star HD4308. Interesting characteristics of this planet 
include its Uranus-like mass and the low metallicity of its par- 
ent star. The paper is structured as follows. Section 2 briefly 
describes the physical properties of the parent star The radial- 
velocity measurements, observation strategy, and orbital solu- 
tion are presented in Sects. 3 and 4. Finally, we discuss the 
characteristics of the new planet in Sect. 5 in the context of the 
latest core-accretion planet-formation models, whereas the last 
section is devoted to some concluding remarks. 

2. Stellar characteristics of HD4308 

The basic photometric (G5 V, V = 6.54, B-V^ 0.65) and as- 
trometric = 45.76 mas) prop erties of H P 4308 were taken 
from the Hipparcos catalogue Tesa" 1997). They are recalled 
in Tabled together with infeiTed quantities like the absolute 
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Fig. 2. CanH (/I = 3968.47 A) absorption Hne region of the 
summed HARPS spectra for HD 4308. Re-emission in the core 
of the line is absent, showing the low chromospheric activity of 
the star. 



magnitude (My - 4.85) and the stellar physical characteris- 
tics derived from the HARPS spectra followiri g the method 
described in "Santo s et alJ ll200li l2004al l2005l) . A standard 
local thermodynamical equilibrium (UTE) analysis was ap- 
plied to a resulting spectrum built as the sum of the 205 
individual spectra gathered for the star (see next section). 
On this resulting spectrum, we measured an equivalent S/N 
of ~ 1100 in the Li wavelength region (6700 A). This study 
thus provides very precise values for the effective tempera- 
ture (Feff = 5686 + 13), metaUicity ([Fe/H] = -0.31 + 0.01), and 
surface gravity (log g - 4.49 + 0.07) of the star. The quoted un- 
certainties do not include systematic errors, such as the use of 
different temperature scales. However, these systematic errors 
should be small, in parti cular the ones conce rning stellar metal- 
licity (see discussion in ISantos et alJ200 4al 2005). 

Very interestingly, HD 4308 has a sub-solar metallicity con- 
trary to most planet-host stars. The high values of the spa- 
tial velocities of the star (f/= -52kms"'; V=-110kms"'; 
W - -29kms"') further points towards the star probably be- 
longing to the thick disk. 

From the colour index, the derived effective temperature, 
and the corresponding bolometric correction, we estimated the 
star luminosity as 0.99 L©. Except for its metallicity, the star 
is thus very similar to the Sun. Because of the low metallicity, 
however, we interpolated a sub-solar mass (M* = 0.83 M©) in 
the grid of Genev a steUar evolutionary m odels with appropriate 
metal abundance (Schae rer et aljfl993l) . The age derived from 
the model also points towards an old star (ago lOGyr). The 
projected rotational velocity vsin; = 1.2kms ' was estimated 
using the calibration of the CORALIE cross-correlation func- 
tion given in .Santos et al. (.2002) . Table [2 gathers those values 
as well. 

We also computed the \ogR'^^ activity indicator from the 
spectra, measuring the re-emission flux in the CaiiH and K 
lines corrected for the photospheric flux contribution. This in- 
dex represents a useful tool for estimating the stellar radial- 
velocity jitter expected for the star due to rotational modula- 
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Table 1. Observed and inferred parameters of HD4308. 
Photometric and astro metric para meters were taken from the 
Hipparcos catalogue ( ESAlll997h . The other stellar physical 
quantities were obtained from a high-resolution ETL spectral 
analysis or were interpolated in the grid of Geneva evolution 
models (ISchaereretalJl993l) . 



Parameter 




HD 4308 


Sp 




G5V 


V 


[mag] 


6.54 


B - V 


[mag] 


0.65 


n 


[mas] 


45.76 


Mv 


[mag] 


4.85 




[K] 


5685 + 13 


^ogg 


[cgs] 


4.49 ± 0.07 


[Fe/H] 


[dex] 


-0.31 ± 0.01 


L 


[Le] 


0.99 


M. 


[Mo] 


0.83 ± 0.01 


V sin ; 


[km s-'] 


1.2 






-4.93 




[days] 


24 


age(model/logi?^^) 


[Gyr] 


> 10/4.3 



tion of star spots or other active regi ons on the stellar surface 
jSaar & Donahuelll997trSantos et al.l (2000). From this indica- 
tor we derived a stellar rotation period Proi - 24 days (follow- 
ing iNoves et al.lll984l) . as well as an activity-based age es- 
timate of ~ 4.3 Gy r (from the calibration by Donahue 1993; 
lHenrvetalJ ll996»)'. HD4308 was found to be non-active as 
shown by the low value of log/?[^j^ - -4.93^ and the absence of 
re-emission observed in the core of the Call H line displayed 
in Fig.|5] Together with the measured small v sin /, these fea- 
tures indicate very low activity-induced radial- velocity jitter for 
this star. IVIoreover, analysis of the bisector shape of the cross- 
corr elation function (see details on the method in lOueloz et alJ 
I2OOI ) shows no variations in the CCF profile down to the pho- 
ton noise level (Fig.|3}, giving strong support to a planetary ori- 
gin of the radial-velocity variations. 

3. Very high precision radial-velocity 
measurements and observing strategy 

HD4308 is part of the HARPS "high-precision" GTO subsam- 
ple that aims at detecting very low-mass extra-solar planets by 
pushing the radial velocity measurement accuracy below the 
1 ms"' mark. However, asteroseismology observations carried 
out by HARPS have made clear that the precision the instru- 
ment is capable of achieving is no longer set by instrumental 
characteristics but rather b y the stars themselves (Mayor et al. 
l2Q03t iBouchv et alJl2005l) . Indeed, stellar p-mode oscillations 
on short time-scales and stellar jitter on longer time-scales can 
and do induce significant radial-velocity changes at the level of 



' Although the rel ation has been shown n ot to be reliable for ages 
greater than ~ 2 Gyr jPace & PasQuinilf2004) . it is still helpful to dis- 
tinguish between active, young stars and chromospherically quiet, old 
stars. 

^ A value logiJ^^ = -5.05 is given in lHenrv et aljfl993) leading to 
an age estimate of 6.3 Gyr. 




Radial velocity [m/s] 

Fig. 3. Bisector velocity span as a function of the measured ra- 
dial velocity. No correlation is found supporting a non-intrinsic 
origin for the radial-velocity variation. 

accuracy of HARPS measurements. For instance, even a very 
"quiet" G or K dwarf shows oscillation modes of several tens 
of cms ' each, which might a dd up to radial-velo city ampli- 
tudes as large as several ms"' jBouchv et al.''2005'). As a con- 
sequence, any exposure with a shorter integration time than the 
oscillation period of the star, or even than mode-interference 
variation time-scales, might fall arbitrarily on a peak or on a 
valley of these mode interferences and thus introduce addi- 
tional radial-velocity "noise". This phenomenon could, there- 
fore, seriously compromise the ability to detect very low-mass 
planets around solar-type stars by means of the radial-velocity 
technique. 

To minimize these effects as much as possible, the stars 
of the HARPS GTO "high-precision" sample ha ve been se- 
lected from the CORALIE planet-search sample dUdrv et alJ 
bOOOl) as slowly rotating, non-evolved, and low-activity stars 
with no obvious radial-velocity variations at the CORALIE 
precision level (typically 10 ms ' for HD4308). Moreover, in 
order to average out stellar oscillations, the observations are 
designed to last at least 15 minutes on the target, splitting them 
in several exposures, when required, to avoid CCD saturation. 
The radial velocities of individual exposures are obtained with 
the standard HARPS reduction pipeline, based on the cross- 
correlation with an appropriate stellar template, the precise 
nightly wavelength calibration with ThAr spectra and the track- 
ing of instrumental drifts w ith the simultaneous ThAr tech- 
nique jBaranne et aljfl99^ . The final radial velocity is then 
obtained by averaging the multiple consecutive exposures. 

This strategy is now applied to all stars in the "high- 
precision" programme. The results are best summarised by the 
histogram of the radial-velocity dispersion of the HARPS mea- 
surements for this programme (Fig.|^. The distribution mode 
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Table 2. Multiple-exposure averaged radial-velocities and un- 
certainties (including calibration errors) for HD4308. All data 
are relative to the solar system barycentre. 



<7rv [m/s] 

Fig. 4. Histogram of the observed radial-velocity dispersion 
((Trv) of the stars in the HARPS "high-precision" sub- 
programme. It has been produced using the data obtained with 
the new observing strategy and includes only stars with more 
than 3 measurements. The position of the planets detected with 
HARPS is indicated by the hatched area. 



is just below 2 ms and the peak decreases rapidly towards 
higher values. More than 80 % of the stars show smaller disper- 
sion than 5 ms"', and more than 35 % have dispersions below 
2 ms ' . It must be noted that the computed dispersion includes 
photon-noise error, wavelength-calibration error, stellar oscil- 
lations and jitter, and, in particular, it is "polluted" by known 
extrasolar planets (hatched part in Fig. 13 and still undetected 
planetary companions. 



4. Orbital Keplerian solution for HD 4308 b 

A set of 205 individual high S/N spectra covering 680 days 
were gathered for HD 4308, as described in Sect|3 They finally 
corresponded to 41 radial-velocity measurements after the 
multiple-exposure averages. Typical exposure times of 5 min- 
utes yielded S/N ratios of ~ 150 per pixel at 550 nm for indi- 
vidual spectra, corresponding to individual photon-noise errors 
of ~ 0.5 ms"'. After combining the spectra, we quadratically 
added the night wavelength calibration error^. The list of the 
averaged radial-velocity measurements obtained for HD 4308 
is given in Table|2l 

The intermediate-season velocities are displayed in Fig.|5l 
with the best Keplerian solution fitted to the data. The derived 



Presently on the order of 80 cms"', this error represents the main 
limitation on the measurement precision, but it will improve soon 
thanks to an ongoing redefiniti on of the thorium lines used for the 
wavelength calibration. See also Lovis et al. 1 2005) for further details 
on the difl^'erent error sources of the radial- velocity estimate 
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period is 15.56 days, and the radial-velocity semi-amplitude K 
is 4.1 ms"'. Despite the small amplitude of the radial-velocity 
variation, we estimated a small false-alarm probability for the 
period derived of 10"^ from Monte-Carlo simulations. The best 
Keplerian planetary solution is circular. Considering a primary 
mass of 0.83 Mq, these parameters lead to a minimum mass 
m2 sin / = 0.0442 Mjup - 14. 1 M^ and a separation a = 0. 115 AU 
for the planet. Figure|6l shows the radial velocities folded to 
the orbital period, including 3 older "single-shot" measure- 
ments obtained before the change in our observational strategy. 
Orbital and infeiTed parameters for the non-circular solution 
are given in Table|3] 
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Fig. 5. Intermediate season of HARPS radial velocities for 
HD 4308. The best Keplerian fit to the data (solid curve) gives a 
minimum mass of 14.1 and an orbital period of 15.56 days 
for the planet. 



Fig. 6. Phase-folded radial velocities of HD 4308 superimposed 
on the best Keplerian orbital solution (solid curve). The 3 older 
points obtained before introducing the new observing strategy 
(see text) are indicated by open symbols. 



Table 3. Orbital and physical parameters for HD 4308 b 



Parameter 




HD 4308 b 


P 


[days] 


15.56 ±0.02 


T 

e 


[JD-2400000] 


53314.70 + 2.0 
0.00 + 0.01 


V 


[km s-'] 


95.2457 ± 0.0002 


a) 


[deg] 


359 ± 47 


K 


[m s-'] 


4.07 ± 0.2 


a[ sin ( 


[10-« AU] 


5.826 


fim) 


[10-"Mo] 


1.09 


1112 sin i 


[Mj„p] 


0.0442 


a 


[AU] 


0.115 


^meas 




41 


Span 


[days] 


680 


cr (O-C) 


[ms-'] 


1.3 


^rcd 




1.92 



The weighted rms of the residuals around the solution is 
1.3 ms which is larger than the internal errors; but, as in 
the case of yu Ara (Santos et al. 2004b), they are most prob- 
ably due to r esiduals of unaveraged stellar oscillation modes 
teouchv et al.k2005.) . No additional trend of the radial veloci- 
ties is expected, as CORALIE measurements show that the star 
is constant at a 10ms"' level over close to 5 years. 

5. Discussion 

5.1. Constraining formation and evolution scenarios 

Low-mass giant planets, i.e. planets with masses in the range 
10- 100 Me, are particularly interesting as they provide poten- 



tially strong constraints on current models of giant planet for- 
mation and evolution. Indeed, and perhaps contrary to intu- 
ition, the formation of these objects within the current the- 
oretical models appears more difficult than the formation of 
their more massive counterparts. For this reason, objects with 
masses within or at the edge of this range, like juArac and 
HD 4308 b, are especially interesting. 

In the direct collapse scenario, planets form through gravi- 
tational collapse of patches of the proto-planetary disk ( Boss! 

on very sho rt timescales. High-resolution simulations 
(lMaveretalJ2002l E004) of this process show that planets tend 
to form on elliptical orbits with a semi-major axis of several as- 
tronomical units and masses between 1 and 7Mjup. In this sce- 
nario, /i Ara-type planets would have to result from the subse- 
quent evolution involving migration and from a very significant 
amount of mass loss. 

In the framew ork of the core accretion model (e.g. 
IPollack etal."l99^, the final mass of a planet is actually de- 
termined by the amount of gas the core accretes after it has 
reached a critical mass, which is about 10-15 Me. In llda & LinI 
(2004a), this amount is determined by the rate at which gas can 
be accreted (essentially the Kelvin-Helmholtz timescale) and 
by the total amount of gas available within the planet's gravita- 
tional reach. Since for super-critical cores the Kelvin-Helmoltz 
timescale is short and the amount of gas available large (com- 
pared to an Earth's mass) even in low mass disks, planets tend 
to form either less or significantly more massive than the crit- 
ical mass. From a large number of formation model calcula- 
tions, [w^^Xi^llloO^ have found that only very few planets 
form in the mass range 10- 100 Me, a range they actually call a 
planetary desert. 
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In the extended core accretion models by lAlibert et alJ 



(|2005|), the planet should be able to accrete gas over the entire 
lifetime of the disk, due to its migration. However, since the 
latter thins out with time and the planet eventually opens a gap 
as it grows more massive, the gas supply decreases with time. 
The growth rate of the planet is actually set by the rate at which 
the disk can supply the gas rather than by the rate at which the 
planet can accrete it. Monte Carlo simulations are ongoing to 
check whether these r nodels lead to a di fferent planetary initial 
mass function than in llda & Linl(l2Q04ah . 



At first glance, the relatively numerous small-mass objects 
discovered so far s eem to pose a pro blem for current planet 
formation theories (^L ovis etaPbOOSb : however, the situation 
is actually more complex. For example, ex amining a poten- 
tial m etallicity-mass relation for exoplanets, iRice & Armitag el 
1I2OO 5) raise the possibility that the planet desert might be popu- 
lated by late-forming planets for which the evolution/migration 
is stopped at intermediate masses and distances by the disso- 
lution of the disk. Also, since all the very low-mass planets 
that are cuiTently known are located close to their star, one 
cannot exclude the possibility that these objects were initially 
much more massive and then have lost a significant amount 
of their mass through ev aporation during their lifetime (see 
iBaraff'e et al .12004 120Q5L for a more detailed discussion). 



While mass loss from initially more massive objects could 
possibly account for the light planets very close to their star, 
it is not clear whether ju Ara c located at a distance of 0.09 AU 
could actually result from the evaporation of a more massive 
object. The situation is even more critical for HD 4308 b that is 
located farther away (0. 115 AU) from its parent star, which is, 
in addition, less luminous than yu Ara by a factor of ~ 1.8. The 
effects could possibly be compensated for, at least partially, by 
the estimated very old age of the star. However, as more fi Ara- 
like objects are being discovered, and if they are all the results 
of the evaporation of larger mass planets, the question of the 
probability of detecting these systems shortly before complete 
evaporation will become a central one. As already stated by 
^araffe et al. ( 2005), the current evaporation models are still 
affected by large uncertainties - lack of detailed chemistry, 
non-standard chemical composition in the envelope, effect or 
rocky/icy cores, etc. - that will need to be clarified in order 
to solve the question of the possible formation of fi Ara-type 
planets through evaporation. 

Finally, given their close location to their star, the detected 
small mass planets are likely to have migrated to their current 
position from further out in the nebula. The chemical com- 
position of these planets will depend upon the extent of their 
migration and the thermal history of the nebula and, hence, 
on the composition of the planetesimals along the accretion 
path of the planet. The situation is complicated by the fact 
that the ice-line itself is moving as the nebula evolves (see e.g. 
LSasselov & L ecar"2000V Detailed models of planetary forma- 
tion that include these effects have yet to be developed. 



5.2. Influence of parent star metallicity 
It is well-established that the detected 



preferentially orbit metal-rich stars 


(Gonzalez 


1998; 


Gonzalez & Vanture 1998; Gonzalez et al. 


2001;lSantos etal] 


.200L .2003. .2005; .Fischer & Valentii i2005). The frequency 



of planets is even found to be a steeply rising function of 
the parent star's metallicity, as soon as the l atter is over solar 
dSantos et alJl200lt Ipischer & Valentii l2005l) . In the scenario 
where gas giants acquire their mass through planetesimal 
coagulation followed by rapid gas accretion onto the core, the 
high probability of a planet to be present around metal-rich 
stars arises naturally"* if protostellar disks attain the s ame frac- 
tion of heavy material as the forming central star dlda & LinI 
l2004bh . 

If the new hot Neptune planets are the remains of evap- 
orated ancient giant planets, they should also follow the 
metallicity trend of their giant progenitors. This does not 
seem to be the case, considering tha t the 7 known plan - 
ets with m?sini<21My - / zArac (Santos etal. 2004^), 

2004 



55Cnc jMcArthur et 



G1777AC 



cArtnureta 
dVogf et alJ 



G1436 (Butler etal. 



2005 



120 05). G1876d (Rivera etal] , 

Gl 58 1 jBonfils et alj2005i) . and HD 4308 b - have metallicities 
of 0.33, 0.35, 0.02, 0.14, -0.03, -0.25, and -0.31, respectively. 
Although the statistics are still poor, the spread of these values 
over the nearly full range of planet-host metallicities (Fig.0 
suggests a different relation between metal content and planet 
frequency for the icy/rocky planets in regard to the giant ones. 

However, we have to note here that 3 of the candidates 
orbit M-dwarf primaries. Recent Monte-Carlo simulations by 
(ida & Lin ( 2005|) show that planet formation around small- 
mass primaries tends to form planets with lower masses in the 
Uranus/Neptune domain. A similar result that favours lower- 
mass planets is also observed for solar-type stars in the case of 
the low metallicity of the protostellar nebula dlda & Linl2004bl 
Mordasini et al. in prep). Future improvements in the planet- 
formation models and new detections of very-low mass planets 
will help to better understand these 2 converging effects. 

6. Summary and concluding remarks 

In this paper we have reported the detection of a new very light 
planet in the Uranus/Neptune mass range that orbits the defi- 
cient star HD 4308. With a period of 15.56 days (a ^ 0.1 15 AU), 
the planet is very similar to /i Arac jSantos et alJl2004'hl) . This 
new discovery was made possible by the high radial-velocity 
precision reached with the HARPS spectrograph and the ap- 
plied observing strategy of averaging down the intrinsic radial- 
velocity "noise" induced by stellar oscillation modes. The 
HARPS precision is demonstrated by the very low residuals 
around the orbital fit, allowing for accurate determination of 
the orbital parameters. 

The existence of small planets (lO-lOOM®) at short and 
intermediate distances from their stars strongly constrains the 
standard planet formation and evolution theories. The difficulty 



On the contrary, such a correlation is not expected with the 
gravitational-instability scenario. Recent simulations even lead to the 
opposite result (.Cai et al..2005i') . 
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-0.2 0.2 

[Fe/H] (all sources) 

Fig. 7. Metallicity distribution of the sample of extrasolar 
planet hosts for planets with shorter periods than 20 days. Stars 
with Neptune-mass planets are indicated by the shaded his- 
togram. 



arises from the fast runaway accretion and the potentially large 
amount of gas available for accretion, both leading towards 
larger mass planets whose inward migration should turn them 
into hot Jupiters (or hot Saturns). However, other effects, such 
as metallicity (and HD4308 is metal deficient), have not yet 
been studied in detail. 

Evaporation could possibly be invoke d to account for the 
plane ts that are very close to their star dBaraffe et al.ll2004l 
I2OO5I) . However, it is not clear whether /i Arab, and especially 
HD 4308 b, could actually result from the evaporation of a more 
massive object. In the case of negligible evaporation or if other 
small mass planets should be discovered for which evaporation 
can safely be ruled out, the existence of the "planetary desert", 
or at least its depth, must be questioned. Therefore, from an ob- 
servational point of view, a larger sample of small objects far 
enough away from the star, would be of paramount importance 
for further constraining our understanding of the formation of 
giant planets. 

Recent HARPS discoveries indicate that a population of 
Neptune- and Saturn-mass planets remains to be discovered be- 
low 1 AU. The increasing precision of the radial-velocity sur- 
veys will help answer this question in the near future, thereby 
providing us with useful new constraints on planet formation 
theories. With the precision level achieved by HARPS, a new 
field in the search for extrasolar planets is now open, allow- 
ing the detection of companions of a few Earth masses around 
solar-type stars. Very low-mass planets (<10Me) might be 
more frequent than the previously found giant worlds. Such 
planets will furthermore be preferential targets for space mis- 
sions like the photometric satellites COROT and Kepler. 
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